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ABSTRACT 

Sub-arcsecond lensing statistics depend sensitively on the inner mass profiles 
of low-mass objects and the faint-end slopes of the Schechter Inminosity function 
and the Press-Schechter mass function. By requiring the luminosity and mass 
functions to give consistent predictions for the distribution of image separation 
below 1”, we show that dark matter halos with masses below 10^‘^Mq cannot have 
a single type of prohle, be it the singular isothermal sphere (SIS) or the shallower 
“universal” dark matter prohle. Instead, consistent results are achieved if we 
allow a fraction of the halos at a given mass to be luminous with the SIS prohle, 
and the rest to be dark with an inner logarithmic slope shallower than —1.5 
to compensate for the steeper faint-end slope of the mass function compared 
with the luminosity function. We quantify how rapidly the SIS fraction must 
decrease with decreasing halo mass, thereby providing a statistical measure for 
the ehectiveness of feedback processes on the baryon content in low-mass halos. 

Subject headings: galaxies: evolution - galaxies: structure - gravitational lensing 


1. Introduction 

The distribution of the luminosity of galaxies and the distribution of the mass of dark 
matter halos are well approximated by the Schechter luminosity function (Schechter 1976) 
and the Press-Schechter mass function (Press & Schechter 1974), respectively. Both functions 
increase as a power law towards the low luminosity and mass ends, but the mass function 
increases with a steeper slope than the luminosity function. Low mass halos must therefore 
contain relatively less luminous baryonic material in comparison with massive halos. Detailed 
models of galaxy formation have been able to account for this difference by feedback processes 
such as supernova explosions, stellar winds, and photo-ionizations that suppress the amount 
of baryons and star formation rates in low mass halos (e.g., Benson et al. 2002; Somerville 
& Primack 1999; Kaufma nn et al 1993). 
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In this Letter we examine this issue from a different perspective of small-separation 
strong gravitational lensing. The image separation distribution of lenses below 1” depends 
sensitively on both the inner mass prohle of galactic halos and the faint end slope of the mass 
and luminosity functions. We compare the traditional approach that models the lenses as SIS 
and the Schechter luminosity function, with a dark-matter based approach that models the 
lenses with a certain halo mass prohle and the Press-Schechter mass function. We investigate 
the constraints on the inner total mass prohles of halos by requiring the two approaches to 
give consistent predictions. Since evidence based on stellar dynamics of elliptical galaxies 
(e.g., Rix et al. 1997; Romanowsky & Kochanek 1999; Treu & Koopmans 2002), modeling 
of lensed systems (e.g., Cohn et al. 2001), and hux ratios of multiple images (Rusin & Ma 
2001; Rusin 2002) all give an inner prohle for lensing galaxies that is consistent with SIS, 
we will use SIS in the lensing calculation with the luminosity function. Dark matter halos 
then clearly cannot all be SIS because if so, the mass function is steeper than the luminosity 
function and would lead to a relatively higher lensing rate at smaller angular scale. We 
will show that modifying the SIS to any single hatter prohle for halos does not work either. 
Instead, we discuss in § 4 how a combination of prohles is needed to resolve the problem. 

This paper complements several recent studies on strong lensing statistics in which the 
emphasis is on the ehects of lens mass prohles and baryon compression on the cumulative 
lensing rates at 6^ > 1" and the implications for cosmological parameters from the scarcity of 
large separation (> 3") systems (Keeton & Madau 2001; Keeton 2001; Kochanek & White 
2001; Kochanek 2001; Sarbu, Rusin, & Ma 2001; Li & Ostriker 2002; Oguri 2002). It is 
pointed out that modeling cluster-scale lenses with mass prohles shallower than the SIS 
greatly reduces the lensing rate and brings the concordance CDM model predictions into 
agreement with observations. The focus here is on the less explored sub-arcsec range. We 
use the predicted shape for the diherential distribution of image separation to quantify how 
rapidly the fraction of SIS halos must decrease with decreasing mass. 

In this paper the cosmological model is taken to have a present-day matter density Dm = 
0.3 (with 0.05 in baryons), cosmological constant Da = 0.7, Hubble parameter h = 0.75, and 
matter huctuation ag = 0.92. The lens potentials are assumed to be spherically symmetric 
because we are mainly concerned with the lensing optical depth, which is more sensitive to 
the velocity dispersion and the radial prohle of the lens than its ellipticity (Kochanek & 
Blandford 1987). The luminosity function is assumed to have a constant comoving galaxy 
number density, which is consistent with the nearly constant comoving halo number density 
(for a hxed velocity) up to redshift ~ 5 in the Press-Schechter formula (Bullock et al. 2001). 
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2. Lensing Rates from Luminosity and Velocity Functions 


The galaxy luminosity function takes the form (Schechter 1976) 

,l,(L)dL = ,j,. (£)”«■'■''■•<'(£) ■ (1) 

An alternative measure is given by the (circular) velocity function 


which is related to the luminosity function by L oc u”, P + 1 = n{a + 1), and 
The velocity function can be derived from galaxy survey luminosity functions and kinematic 
luminosity-velocity relations, e.g., from various large pre-SDSS optical surveys (Gonzalez et 
ah 2000) and the 2MASS infrared survey (Kochanek 2001). The image separation distribu¬ 
tion of lenses at angle 6 is related to the galaxy velocity function by 


dP 

le 


Zl) aiens{Vc, Zl) B , 

dzi dO 


(3) 


where dr/dz = cHp^{l + + zp + (for 12^ + a = 1), o'lens is the lensing 

cross section, B is the magnihcation bias, and iplvc, zpdvc is the physical number density of 
galaxies with circular velocity between Vc and Vc + dvc at lens redshift zi. 

An SIS lens has a density prohle p(r) = v‘^/{8nGr‘^) and produces an image separation 
of 6* = 29e where 9e = ‘2tt{vc/cPDis/D g is the Einstein radius. For 9 = 1" and source 
redshift Zg = 1.2, Vc ranges from 225 to 325 km/s for zi = 0.3 to 0.7. The lensing cross 
section is uiens = t^{9eDiY = At[^{vc/ cY{DiDis/Dgp (Schneider, Ehlers & Falco 1992). For 
small 9 (< 1"), one can show analytically that the slope of dP/d9 in eq. (3) depends only on 
the faint-end slope P of the velocity function: 


^ OC ^ ^n(a+l)/2+l ^ 


for small 9 , 


(4) 


where P + 1 = n{a + 1) is used to relate P to the faint-end slope a of the luminosity function 
and the luminosity-velocity relation L oc . The distribution dP/d9 therefore has a positive 
slope on sub-arcsec scale if the velocity function is shallower than P = —3, or if the luminosity 
function is shallower than a = —{1 + 2/n). Fig. 1 illustrates this dependence. Recent surveys 
favor P ~ —1.3 to —1.0 (e.g., Gonzalez et al. 2000; Kochanek 2001), indicating a positive 
slope for dP/d9 at small 9. The image separation distribution of 13 lenses found in the 8958 
GLASS radio sources (Browne et al. 2002) is shown in Fig. 1 for comparison. Note that the 
smallest and largest 9 bins contain only one lens each. 
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3. Lensing Rates from Mass Functions 

Lensing probes mass after all, so let us take an alternative approach by modeling the 
lenses as a population of dark matter halos with a Press-Schechter type of mass fnnction. 
Similar to eq. (3), the image separation distribution is given by 

dP f , dr dM ,^ \ 

le^j ^lens(M, zi) B , (5) 

where n{M, zi)dM is the physical number density of dark halos with mass between M and 
M + dM at zi- We use the improved version of the mass fnnction by Jenkins et al. (2001). 
Unlike the lensing galaxies in eq. (3), the lensing halos in eq. (5) may or may not host 
central baryons depending on if the lens is baryon or dark matter dominated near its center. 
We will therefore consider different inner mass prohles. For luminous lenses in which baryon 
dissipation controls the inner density, we use the SIS prohle as in § 2. For dark lenses withont 
a significant amount of baryons, we consider the shallower prohles fonnd in high resolution 
dark matter simulations: p oc r~^ (Navarro et al. 1997) and r~^'^ (Moore et al. 1999). 

The lensing properties of the three prohles - SIS, NFW, and Moore - are as follows. 
For SIS, we relate the circular velocity Vc of galaxies to the virial velocity Vmr of dark 
halos by Vc = 'jvVvir, where ~ 1.3 to 1.8 from various baryon compression models and 
observational constraints (Oguri 2002 and references therein). We take Oguri’s best ht 
value jjj = 1.67 here. Using v^ir = (47rG^M^pA^jr/3)^'^®, we then have 9 = where 

A = 47rG(jy/c)^(Dis/Ds)(4:7rpAyir/3)^^^ and A^ir zz 178. 

For the Moore et al. prohle p(x) = pS/{x^^‘^ + x^), we hnd the projected snrface 
density S well approximated by k{x) = S/Sc = 3Akq/[x^^‘^{1 + 1.8x)^/^] where x = r/r^. 
Here k is the convergence. Sc = /4tiG){Ds/D iDis) is the critical snrface density, and 
^0 = rspS/T^c- The scale radius is related to the concentration parameter by c = r^ir/vs, 
where r^ir is the halo virial radius, and 5 = 100c^/ln(l + c^/^). The rednced dehection 
angle, related to k by a{x) = 2x~^ dy yK,{y), then has the simple analytic form a{x) = 
(12Ko/v^a;)[ln(v^ + Vl + bx) — + bx], where b = 1.8. The lensing cross section 

is (Tiens = 'n'{Prad.Di)'^, where Prad U the angnlar size of the radial caustic, which we obtain by 
solving the lens equation. 

For the Moore prohle, we hnd the htting fnnction (PradDi) = Q-3rsK,l/{l — I.Ikq^ + 
4.5kq®) accnrate (with < 5% error for kq < 6) and usefnl in speeding np the computation. 
Since the angnlar separation of the outermost images is insensitive to the location of the 
sonrce (Schneider et al. 1992), we use the size of the tangential critical curve for the image 
separation: 6 = 26tan- We hnd the htting fnnction (OtanDi) = 38rsKl/{l — I.Skq® + 19kq^^) 
accnrate (with < 4% error for kq ^ 6). Similarly for the NFW lenses, we use (PradDi) = 
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-4r^fi;o(l - 0.4K[j-i + 0.5k[]-2)/ exp[(3 + Kq ^)/2] for aiens = T^iPradDif, and {OtanDi) = 2r^(l - 
0.7 kq^ + 1.35kq®®)/ exp[(l + Kq^)/2] for 9 = 29tan- 

We compute the magnification bias B from the htting formula given by eq. (21) in Oguri 
et ah (2002). We have tested this formula against numerical calculations and found good 
agreement. A similar £t given by eq. (67) in Li & Ostriker (2002), however, substantially 
underestimates the bias for kq ^ 1 and overestimates it for kq > 1. This is because their fit 
assumed da/dx = 0 and therefore neglected a factor containing (1 — da/dx), where a is the 
deflection angle. We find this not to be a valid assumption in general. 

Fig. 2 compares dP/dd for SIS, Moore, and NFW lenses computed from the halo mass 
function in eq. (5), and dP/dO for SIS lenses computed with velocity functions of different 
slope (3 in eq. (3). It shows that no single mass profile with the halo mass function can 
match the dP/d9 predicted by the observed velocity function of ~ —1.3. The SIS (solid) 
and Moore (long-dashed) profiles predict wrong shapes for dP/d6, a reflection of the steeper 
faint end of the mass function compared with the luminosity function. The shape of dP/dO 
for the shallower NFW prohle (short-dashed) resembles more closely the velocity function 
prediction, but the lensing amplitude is miniscule. We note that the magnification bias B 
has been included in Fig. 2, which is generally significantly higher for shallower inner mass 
profiles, but the resulting NFW lensing amplitude is still much too low. 


4. Mass vs. Light: Resolution 

To bring the predicted shape for the image separation distribution from eq. (5) into 
agreement with eq. (3), we explore the possibility that at a given mass, a fraction of the 
lenses is luminous, baryon dominated at the center and has the SIS prohle, while the rest of 
the lenses is dark matter dominated and has a shallower inner profile. Instead of this bimodal 
model, one can presumably allow the slope of the inner mass profile to decrease smoothly 
with mass. Observations of lensing galaxies, however, consistently find that the combined 
stellar and dark matter mass profile inside the Einstein radius is well fit by the SIS (e.g., Rix 
et al. 1997; Romanowsky & Kochanek 1999; Cohn et ah 2001; Treu & Koopmans 2002). At 
the same time, galaxy formation models show that the energetics of feedback processes are 
sufficient to expel baryons in some < IO^^Mq halos. The bimodal model combining SIS and 
dark matter prohles therefore appears physically motivated and will be used below. 

From the solid and long-dashed curves in Fig. 2, we conclude that a combination of SIS 
and Moore profiles cannot reproduce the shape of dP/dO predicted by the velocity function 
with (3 ~ —1.3. This is because the inner slope of the Moore profile is close enough to SIS 
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that the two predict similar shapes for dP/dd. Making some lenses dark with r~^'^ inner 
prohle will therefore not reproduce the monotonically rising dP/dd at small 9 for the velocity 
function. If halos have the shallower r~^ prohle, however, the dark lenses will have negligible 
lensing optical depth compared with the SIS (short-dashed vs. solid in Fig. 2). We can then 
match the two predictions by parameterizing the fraction of SIS halos at a given mass with 


fsis{M) oc 


^{M/Mc)^2 _ I ’ 


( 6 ) 


which grows as a power law, fsis ~ at small M and falls exponentially at large 

M. This makes the halos mostly SIS on galactic mass scale Me where baryon dissipation is 
important, and mostly NFW on cluster and sub-galactic scales where dark matter dominates 
the potential. Our main interest here is in determining the slope r]i, which has the convenient 
property that it depends only on the faint end slope (3 and not on other parameters in the 
velocity function in eq. (2). It also gives a simple parameterization of the importance of 
feedback processes on the density profile as a function of halo mass. We note that since 
the relation between the image separation 9 and halo mass M is redshift-dependent, the 
factor fsis must be included inside the integral of eq. (5). We do not consider explicit 
redshift dependence in fsis here, which can be put in at the expense of introducing more 
parameters. Future work combining fsis determined from lensing with galaxy formation 
models and simulations may offer useful constraints on the time evolution of fsis- 

Fig. 3 shows the excellent agreement between the two predictions for the shape of dP/d9 
for four faint-end slopes of the velocity function. The required rji in eq. (6) is ~ 0.85, 0.75, 
0.53, and 0.2 for f3 = —1, —1.3, —2, and —3, respectively. The other two parameters Me and 
rj 2 depend on the shape of the velocity function. For n* = 250 km s“^ and n = 2.5 from the 
SSRS2 sample (Gonzalez et al. 2000), we find a good match with Me ~ 4.5 x IO^^Mq and 
ri 2 ~ 0.72. We also hnd it necessary to lower the overall amplitude of dP/d9 for the dashed 
curves in Fig. 3 by ~ 20% to 40% (for normalized fsis shown in Fig. 4) to match the dotted 
curves. We have not attempted to hne tune it since the amplitude of dP/d9 depends on 
several uncertain parameters, e.g., the source redshift distribution, the normalization and 
redshift evolution of the luminosity function, and the precise value of = Ve/v^ir for SIS 
halos. Instead we have focused on the constraints from the shape of dP/d9. 

Fig. 4 shows the required fsis{M) for each of the four velocity functions in Fig. 3. 
Our fsis at the high mass end for f3 = —1.3 agrees well with the result from Oguri (2002), 
which used fsisiM) = exp[l — {M/M^Y^] for M > M^ following Kochanek (2001), and set 
fsis = 1 for M < Mh, i.e., it ignored the dark lens fraction for small masses. By contrast, 
our form of fsis in eq. (6) is a smooth function and takes into account all masses. 

Fig. 4 illustrates that a very steep faint-end slope {(3 ~ —4) for the velocity function will 
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be required if low-mass halos all have the SIS prohle. Current galaxy surveys favor a much 
shallower faint end slope of ~ —1.3. We thus conclude that the percentage of halos that 
can have the SIS prohle must decrease rapidly with decreasing halo mass below This 

implies that feedback processes are increasingly effective in reducing the baryon content in 
small objects, a trend consistent with semi-analytic galaxy formation models. Moreover, the 
halos that have non-SIS prohles must have an inner density of p ~ r~^ or shallower so as to 
contribute negligible lensing optical depth. (The possibly shallow prohles of dwarf galaxies 
will therefore not ahect the lensing predictions here.) The steeper p ~ would predict a 
shape for the lensing image separation diherent from the observed luminosity function. 

I have enjoyed discussions with Paul Schechter, Joe Silk, Aveshi Dekel, David Rusin, 
Dragan Huterer, and David Hogg. This work is supported in part by an Alfred P. Sloan 
Foundation Fellowship, a Cottrell Scholars Award from the Research Corporation, and NASA 
grant NAG 5-12173. 
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6 (arcsec) 

Fig. 1.— Image separation distribution dP/dO calculated from eq. (3) using the galaxy- 
velocity function in eq. (2) and SIS profile. For small 6, the slope of dP/dO only depends on 
the faint-end slope (3 of the velocity function, changing from positive to negative -when the 
velocity function steepens beyond (3 = —3. Other velocity function parameters used here 
are n* = 250 km s“^, n = 2.5, = 0.073/i^ Mpc“^ from the SSRS2 sample in Gonzalez et 

ah (2000). The sources are put at the mean redshift 1.27 of the CLASS radio sources -with 
a power-la-w flux distribution of slope -2.1. The histogram sho-ws the 13 CLASS lenses and 
the la poisson errors (Bro-wne et ah 2002). 
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6 (arcsec) 


Fig. 2.— Image separation distribution dP/dd calculated from galaxy velocity function 
(eq. (3); dotted) vs. mass function (eq. (5)). The result depends strongly on the choice of 
the function and the assumed lens mass profile. The lensing amplitude decreases rapidly as 
the inner mass prohle is lowered from —2 (SIS; solid), to —1.5 (Moore; long-dashed), to —1 
(NFW raised by a factor of 10 to fit in the plot; short-dashed). Note that no single halo 
prohle can bring the mass function prediction for dP/dO into agreement with the velocity 
function prediction. The histogram shows CLASS data as in Fig. 1. 
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6 (arcsec) 

Fig. 3.— Different predictions from velocity vs. mass function in Fig. 2 can be brought 
into agreement if at a given mass, a fraction fsis{M) (eq. (6)) of the dark matter halos is 
assigned SIS and the rest NFW. The predicted dP/dO from the velocity function (dotted) 
and the mass function (dashed) then agree very well for suitable choices of parameters r]i, 772 
and Me for fsis (see text). The histogram shows CLASS data as in Fig. 1. 








M (Mo) 


Fig. 4.— The fraction fsis{M) of halos with SIS prohles needed for the consistent predictions 
in Fig. 3. As the faint-end slope (3 of the velocity function steepens, a larger fraction of low- 
mass halos is allowed to be SIS. Galaxy surveys favor f3 ~ —1.3 (solid), requiring fsis ~ 
with Tji ^ 0.75 below The dotted curve shows Oguri’s result (2002), which agrees 

well with our solid curve at large mass but assumes fsis = 1 at small mass. 




